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Abstract 

We propose a model to describe the concentration dependence of the 
viscosity of soft particles. We incorporate in a very simple way the softness 
of the particles into expressions originally developed for rigid spheres. 
This is done by introducing a concentration-dependent critical packing, 
which is the packing at which the suspension looses fluidity. The resultant 
expression reproduces with high accuracy the experimental results for 
suspensions of star polymers in good solvents. The model allows to explain 
a weak increase of the viscosity observed in the case of diblock copolymer 
stars suggesting that the reason for this peculiar behavior is mainly a 
consequence of the softness of the particles. In the semi-dilute regime, 
suspensions of star polymers are modeled using the Daoud-Cotton picture 
to complete the description in the whole concentration regime. 

It is of basic and applied importance to understand the flow behavior of flu- 
ids containing colloidal particles and macromolecules. The range of industrial 
applications of such complex fluids range from motor oils, coatings, drilling flu- 
ids and food products among others. Many colloidal suspensions consist of hard 
particles whose rheological properties have been studied profusely. By contrast, 
much less attention has been paid to suspensions containing soft colloids. Rep- 
resentatives of soft colloids include deformable particles and particles that in 
addition may interpenetrate to some extent. 

The interest in studying soft particles stems from the desire to explore the 
behavior of a number of complex fluids composed of nonrigid structures such 
as polymerically stabilized colloidal spheres [T], block copolymer micelles [2 , 
star polymers [3], hard spheres with a grafted polymer brush g], dendriti- 
cally branched polymers [5] and others. In particular, star polymers which 
are branched macromolecules with a small central core from where / arms of 
identical linear polymer chains, each one containing N segments emerge are 
considered an important prototypical model of soft particle because their size 
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and architecture can be tailored between that of linear polymers and ultra-soft 
colloids to nearly hard-spheres [B]. 

Pioneering theoretical work to calculate the viscosity of dilute colloidal sus- 
pensions was initiated by Einstein who obtained an expression to calculate the 
viscosity as function of the concentration of dilute suspensions of spherical hard 
particles [7] . Various models to calculate the viscosity of dilute suspensions of 
particles other than solid spheres have been proposed. Among them we can 
mention the case of emulsions of spherical droplets [5] , homogeneously porous 
rigid spheres [9]- [IT] and uncharged spherical soft particles |12l. 

For non-dilute systems, the suspension rheology is determined by the inter- 
play between the direct particle-particle interactions and the solvent-mediated 
hydrodynamic interactions (His) 1 1 3] . Their many body nature impose a formidable 
difficulty for the calculation of rheological quantities such as the shear viscosity 
of the system. The fact that many particles are soft complicates even more the 
calculation. Nonetheless, it is a common practice to correlate the concentration 
dependence of the viscosity of star polymers and other hairy particles using 
phenomenological models designed for solid hard spheres [T3], [15]. In partic- 
ular, experimental data are frequently compared with the expression proposed 
by Krieger and Dougherty for hard colloidal particles [16_ 



where rj (4>) is the viscosity of the suspension, 770 the viscosity of the background 
solvent, (j> is the volume fraction of the colloidal particles, ^ max is the volume 
fraction at maximum packing and [77] = 5/2 is the intrinsic viscosity of hard 
spheres. Comparing experimental data for star polymers with this expression 
might be accurate for stars with a large number of short arms since they behave 
essentially as hard spheres. However, such comparisons do not give satisfactory 
results in the case of stars with larger arms because in this case the particles 
are much softer. In addition to this, it has been shown that the Krieger and 
Dougherty expression underestimates the viscosity of the suspension at large 
volume fractions even in the case of hard spheres, as explained thoroughly in 



Alternative expressions have been devised to predict the viscosity for concen- 
trated suspensions including in a more accurate way hydrodynamic and excluded 
volume interactions. Using a differential effective medium technique based on 
a progressive addition of particles to the sample in which the new particles in- 
teract in an effective way with those added in previous stages, we have shown 
[17 -[20 that the static and high-frequency viscosities of suspensions of colloidal 
particles can be correlated with the expression 



where [rj] is the intrinsic viscosity and k a fitting constant related to the critical 
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Ref. [TF]. 
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volume fraction <f> c at which the suspension loses its fluidity and is given by 




Expression (|2]) incorporates crowding effects and hydrodynamic interactions 
[17 -[20 and has been applied successfully to suspensions of hard spheres |18| . 
emulsions of spherical droplets [19] , suspensions of arbitrarily-shaped hard parti- 
cles [T7] , suspensions of rigid core-shell permeable particles [5U] and suspensions 
with power-law matrices [51] with excellent results. However, as in the case of 
the Krieger and Dougherty model, when applied to soft (interpenetrable) parti- 
cles, Eq. Q fails to reproduce the viscosity-concentration curves satisfactorily. 

The purpose of the present work consists in propose a method to incorporate 
the softness of the particles in the above mentioned expressions to obtain a model 
that better correlates with experimental results. 

Since we have in mind to apply the model to suspensions of star-branched 
polymers we have to take into account the possibility that the particles are 
permeable to some extent to the surrounding solvent. To do so, we consider a 
core-shell particle consisting of a spherical hard core of radius a coated with a 
surface soft layer of thickness <5. The soft layer may represent a grafted polymer 
brush [55] or the outer sections of a star polymer, as represented in Fig. [I] 
Thus, the polymer-coated particle or star polymer has a non-porous inner radius 
a and an a permeable outer shell of radius b = a + S. The volume fraction <p 
is defined as the volume occupied by the particles including the shell, divided 
by the volume of the sample. To model the permeability of the outer shell, the 
polymer segments are regarded as resistance centers distributed in the permeable 
polymer layer, exerting frictional forces — 711 on the liquid flowing in the layer, 
where u is the liquid velocity relative to the particle and 7 is the frictional 
coefficient. T he res ult of this model is [15] [77] = (5/2)L 2 (Xb,a/b) jL x (Xb,a/b), 
where A = y/j / 770 . Expressions for L\ (Xb,a/b) and L2 (Xb,a/b) can be found 
in Ref. [50] ■ Therefore, the permeability of the outer shell is incorporated in 
the intrinsic viscosity [77] of the core-shell spheres. 

If we consider that the outer core of the soft sphere do not exert frictional 
forces on the liquid flowing through it, that is, A = 0, then the expression for 
[77] simplifies considerably and one gets 

H-IG)'- W 

The present model, Eq. ^ together with Eq. Q allows to obtain the viscos- 
ity as a function of the concentration of rigid porous core-shell particles and has 
been successfully applied to calculate the static viscosity of surfactant coated 
particles and the high-frequency viscosity of homogeneous porous particles (20] . 
However it does not give satisfactory results for suspensions of star polymers. 

As mentioned above, experimental results of the concentration-dependent 
viscosity of star polymers can not be satisfactorily reproduced with models 
developed for hard colloids like the Krieger-Dougherty model, Eq. (fTl or the 
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one given by Eq. (J2J) even when the intrinsic viscosity for core-shell particles, 
Eq. Q, is used. The problem with these expressions is that they assume that 
the particles are rigid. It is apparent that the softness of the stars must be taken 
into account to obtain better results. Here we propose to model the softness of 
the particles by assuming that the critical volume fraction, <p c , is larger than the 
random close packing value 4>rcp since soft particles can interpenetrate each 
other. Specifically, we propose the following ansatz 

4> c = <t>R C p + ^ a , (5) 

where a and /3 are unknown constants that should be determined from the 
fitting to the experimental data. Expression ^ can be interpreted as follows, 
when the volume fraction is very low, the soft particles do not interact too 
much with each other and therefore they do not overlap too much on average. 
In this situation they behave as rigid impenetrable spheres and <j) c — 4>rcp- As 
the volume fraction is increased, the particles are forced to collide more often 
increasing the average overlap which in turn increases the value of <j) c . Eq. (|5| 
then proposes that this increase follows a power law on the volume fraction. 

To test the correctness of this proposal we compare the predictions of our 
model given by Eq. Q including the possibility that the outer sections of the 
star are solvent permeable, using the intrinsic viscosity Q. Finally and most 
importantly, we incorporate the softness of the particles substituting Eq. (|5| 
in Eq. ^ . This is done in Fig. [2] where we compare the predictions of our 
model with the experimental results of Roovers, Ref. [TS] for the relative vis- 
cosity of polybutadiene star polymers in the good solvent toluene. Typically, 
the viscosity vs concentration curves for star polymers are plotted in terms 
of the polymer concentration c/c*, where c* is the overlap concentration de- 
fined by c* = [(N A /M) (4?r/3) , where N A is Avogadro's number, M the 
molecular weight, and Ry the viscosimetric equivalent hard-sphere radius of the 
polymer coil. On the other hand, models for the viscosity of colloidal suspen- 
sions like Eq. {2} are written in terms of the particle volume fraction (f>. In the 
dilute regime these two quantities are related by 

<p = c/c\ (6) 

which allows to plot our model in terms of the polymer concentration. In Fig. 
[2] we perform such plot. Panel (a) shows the experimental results for stars with 
/ = 128 arms. Notice that the viscosity-concentration curve is not very different 
from the hard-sphere result at low concentrations. However at large concentra- 
tions a slight departure from the hard sphere behavior is observed. Our model 
reproduces very well this departure with the fitting parameters shown in Table 
[l] If the number of arms is decreased, one expects a smaller steric repulsion 
between the arms of different stars allowing a larger overlap and therefore a 
larger departure from the hard-sphere behavior. This is shown in Fig. [2] (b) 
where results for a star with 64 arms are shown. The softness of the star is 
larger then in the case of the star with 128 arms and the departure from the 
hard-sphere behavior is larger. Our model predicts correctly the experimental 
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data and the increased softness of this case is reflected in a larger value of the 
parameter /?, as shown in Table [l] Panel (c) shows the results for a star with 
32 arms. The low number of arms of this case allows a much larger overlap 
allowing that concentrations above c* can be reached. Experimental data show 
clearly a different behavior below and above c*. The region below c* can still 
be described with our model with an even larger value of f3 than in the previ- 
ous cases. However, above c* , the overlap is so large that the conformation of 
the star polymers should be described differently. Here we use the Daoud and 
Cotton picture [23] to describe the star conformation in this semi-dilute regime. 

In the Daoud and Cotton model [23], a single star polymer is regarded as 
a spherical region of a semi-dilute polymer solution with a local, position de- 
pendent screening length £ (r), where r is the distance from the center. This 
is represented pictorially by associating with each arm a string of blobs of in- 
creasing size £ (r) as shown in Fig. [3] The blob size varies as £(r) ~ r/ -1 / 2 
and the corresponding local polymer volume fraction as c s (?*) ~ / 2 / 3 (b/r) 4 ^ 3 . 
The size of the star can be obtained from monomer conservation (24] to obtain 
R ~ fttipNf 1 / 5 , where £ p is the Kuhn length. Above the overlapping con- 
centration c* = fN£p I R? , the outer sections of the arms overlap and build a 
semi-dilute solution, often called the sea of blobs, where the inner parts of the 
actual stars are embedded and where the polymer concentration of the sea of 
blobs levels off (see Fig. |3j. The stars embedded in the sea of blobs occupy 
a volume fraction <j> and the sea of blobs the remaining fraction 1 — 0. Above 
the overlapping concentration c* , the values of the polymer concentration in the 
sea of blobs c p , of the fraction occupied by the stars cj>, and of the radius of the 
embedded stars R s d can be calculated from mass conservation |24j : 

/ *\ 5/4 

<P = ki[-) , (7) 
c p ~ c, (8) 



and 



Rsd = Rk 2 ( - J ■ (9) 



c 



The constant k\ can be obtained from continuity of the volume fraction at 
= 4>RCP = 0.637 from Eqs. ([6| and ^ to get k\ = 0.3625, while the constant 
&2 is obtained by continuity from Eqs, Q and Q since R s( j — R at the same 
point, to get k 2 — 0.7130. Notice that an increase of the polymer concentration 
c leads to an increase of the fraction of space occupied by the sea of blobs and 
a shrinkage of the inner star dimensions. Additionally, the stars are embedded 
in a solvent with viscosity T] s b formed by the sea of blobs and from Eqs. Q and 
|9]) the intrinsic viscosity of the embedded stars is 
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The sea of blobs exerts a strong opposition to flow because movement proceeds 
through a disentanglement of the stars from their neighbors [■25] and since they 
consist of a semi-dilute polymer suspension we suggest to model the viscosity 
of the sea of blobs using the scaling relation 

vs b =v° b (~y, (ii) 

where 77^ and 7 are additional fitting constants that allow to model the large 
viscosity of the sea of blobs. 

Finally, substituting Eqs. Q and ( 10 ) and replacing the viscosity of the 



solvent 770 in Eq. ([2| with the viscosity of the sea of blobs rj s b, we obtain the 
viscosity for c above c*. 

Using this approach we can reproduce the viscosity-concentration curve of 
the star with / = 32 above c* as shown in Fig. |2jc). Notice that in this regime, 
the volume fraction occupied by the stars decreases with concentration and 
this effect alone would produce a decrease of the viscosity with concentration. 
However, this decrease is offset by the increase of the viscosity of the sea of 
blobs with concentration. The parameters a and /3 used in this region are the 
ones obtained in the dilute regime and r]® b and 7 are also shown in Table [l] 
Finally, in panel (d) we show the case of a diblock copolymer star with 64 arms, 
an inner core composed of polybutadiene arms with N inner = 20000 and an 
outer shell of polystyrene arms with N outer — 10000. The viscosity data were 
collected in toluene which is a good solvent for both blocks. Again, two different 
regions are apparent and additionally a shoulder appears at a concentration 
near c*. Surprisingly, our model reproduces correctly this shoulder. Therefore, 
our model suggests that the origin of this feature is an intrinsic feature of the 
softness of the particles. The semi-dilute region is again well described by the 
Daoud-Cotton picture. Notice that in all the cases considered the thickness of 
the porous layer is very small as compared to the core size (see the values of the 
parameter a/b in Table [T]). In other words, for any practical purpose, the star 
polymers behave as non porous particles and the value [rj] =5/2 can be used. 

Summarizing, in this work we have presented a simple model for the calcu- 
lation of the concentration dependence of the viscosity for suspensions of soft 
spherical particles consisting of a hard core of radius a covered with a porous 
shell of thickness 5 = b — a. The model contains three main ingredients, the first 
one consists in using an accurate model for the viscosity originally developed for 
rigid particles, Eq. ^ . Secondly, we use a simple expression for the intrinsic vis- 
cosity of permeable core-shell particles, Eq. ([4} to allow the possibility that the 
outer region of the particle is solvent permeable. Third and most importantly, 
we incorporated in a very simple way the effects of the softness of the particles 
through Eq. ^ . Using this model we were able to reproduce the experimental 
results for the viscosity of star-branched polymers and explained the appear- 
ance of a shoulder seen in the viscosity concentration curves of the diblock star 
near the overlap concentration. We conclude that this peculiar behavior is an 
intrinsic feature of the softness of the particles and there is not need to invoke 
alternative mechanisms as a shrinking of the polybutadiene block |15| when the 
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total polymer concentration increases. Additionally, in the semi-dilute regime 
c/c* > 1, the suspension is modeled using the Daoud-Cotton picture which com- 
pletes the description of the viscosity in the whole concentration regime. We 
have tested our model comparing with experimental results finding an excellent 
agreement in spite of the numerous simplifications used in the obtention of our 
model. 
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System 


128 arms 


64 arms 


32 arms 


Diblock star 64 arms 


a/b 


0.988 


0.989 


0.993 


0.984 


a 


1.90 


1.60 


1.74 


1.50 


P 


0.244 


0.580 


0.966 


0.662 


77° 
'Isb 






9.28 


17.4 


7 






2.67 


2.88 



Table 1: Fitting parameters of the model for the systems of Fig. [2] 



Figure Captions 

Fig. 1. Schematics of a model for a suspension of core-shell spheres. A hard 
core of radius a covered by a porous layer of thickness S. The total radius of 
the particle is b = a + 6. 

Fig. 2. Relative viscosity rj{<j))/rjQ as a function of the particle volume 
fraction <f> for star polymers with (a) 128, (b) 64, and (c) 32 arms and (d) diblock 
copolymer micelle with 64 arms in the good solvent toluene. Experimental data 
taken from Ref. |15j . 

Fig. 3. Daoud and Cotton model for a semi-dilute suspension of star poly- 
mers. 
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Figure 1: Schematics of a model for a suspension of core-shell spheres. A hard 
core of radius a covered by a porous layer of thickness 5. The total radius of 
the particle is b = a + 6. 
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Figure 2: Relative viscosity ^(</>)/?7o as a function of the particle volume fraction 
<t> for star polymers with (a) 128, (b) 64, and (c) 32 arms and (d) diblock 
copolymer micelle with 64 arms in the good solvent toluene. Experimental data 
taken from Ref. |15j . 
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Figure 3: Daoud and Cotton model for a semi-dilute suspension of star polymers. 
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